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(57) Abstract 



A multiuser detector for a CDMA system for real-time processing. Interference information of cross correlation from the signals 
of different users is processed to enhance the signal-to-noise ratio by suppressing the interference noise in a mult.ple-stage processmg. A 
sampling processor is implemented to obtain at least two sets of data for each user with different integration umes. The operation of the 
detector is dynamically adjusted to the operating condition of the CDMA system in a near real-time manner. 
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A MULTIUSER DETECTOR FOR CODE DIVISION MULTIPLE ACCESS 
SYSTEMS 

This application claims the benefit of the U. S. 
Provisional Application No. 60/007,469, titled "Multiuser 
5 Detector for Multiple Access Communication System'' , filed on 
November 22, 1995, and 60/007,475, filed November 22, 1995, 
the disclosures of which are incorporated by reference. 

FTELD OF THE INVENTION 

The present invention relates to multiple access 
10 communication systems, and more specifically, to a detec- 
tion system for code division multiple access systems. 

BACKGROUND AND SUMMARY OF THE INVENTION 

Code Division Multiple Access (CDMA) has emerged as a 
key technology for communication services. CDMA 
15 intentionally uses a transmission bandwidth that is 
usually several orders of magnitude larger than the 
minimum required signal or information bandwidth, e.g., 
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the bandwidth implied by the basic message symbol rate in 
a communication system. The power of each user's signal 
is spread over a wide bandwidth. This results in low 
power spectral density and thereby the interference to 
another narrow band signal occupying the same frequency 
range is reduced. It also makes the presence of the 
signal less detectable to an invader. 

A digital signal from each user in a CDMA system is 
modulated with a pseudo-noise (PN) binary sequence that is 
unique to that particular user. This modulation causes 
the spreading over the wide bandwidth. Each PN sequence 
appears random to a naive observer but can be reproduced 
in a deterministic manner by an intended receiver. Any 
two PN sequences are made substantially orthogonal to each 
other but some degree of cross correlation still exists. 
The mutual interference in the same frequency range of 
multiple users is greatly reduced in a CDMA system. This 
orthogonality allows multiple access within the same 
frequency spectrum and makes CDMA systems less vulnerable 
to intentional or unintentional interference. However, 
the cross correlation becomes limiting the system capacity 
as the number of users increases. 

Detection of a CDMA system usually involves cross 
correlation with a locally-generated version of the PN 
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sequence. A desired user signal is generally detected by 
cross correlating with the exact same PN sequence that is 
uniquely assigned to that particular user by the system. 
Signals from other users contribute a small amount of 
5 wideband noise due to non-zero cross correlation value 
produced by the cross correlation operation. 

The ratio of the transmission bandwidth over the 
information bandwidth indicates the degree of the spectrum 
spreading of a CDMA system. This ratio is often referred 

10 as the processing gain. Many beneficial properties of 
CDMA relate to the processing gain. 

FIG. 1 shows a simplified CDMA system 100. User 
messages are modulated with a carrier frequency and 
decoded by PN sequence waveforms prior to transmission. A 

15 receiver system 130 process the received signals by 
filtering and correlation operation to recover the 
embedded user messages. A CDMA detector 137 is an 
essential part of the receiver 130. 

CDMA systems have many unique properties. For 

20 example, multiple access communication by a large number 
of users in a common frequency range in the same and 
neighboring geographical areas are possible. This has 
particular significance in satellite and cellular 
communication systems. The number of simultaneous users 
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is usually proportional to the processing gain of the CDMA 
system. The efficiency of bandwidth usage of CDMA 
increases as the number of simultaneous users grows. The 
resistance to interference increases with the processing 
gain by suppressing the power in a narrow band. Adverse 
effects of multipath signals can be significantly reduced 
due to the spectrum spreading since only a small portion 
of signal will undergo fading at any given time. In 
addition, position location and velocity estimation can be 
measured with high accuracy in a CDMA system. 

The application of CDMA in communication has been 
demonstrated in a number of systems in use including the 
IS-95 based CDMA cellular system disclosed in "On the 
Capacity of A Cellular CDMA System", by Gilhousen et al 
IEEE Trans. Veh. Technol., Vol. VT-40(2), pp303-312, May, 
1991, a PCS system based on IS-95, and the Globalstar LEO 
satellite system. The high capacity, high resolution 
time-of-arrival for positioning measurement, low 
probability of interception, and independent interference 
rejection makes CDMA communication systems one of the most 
promising candidates for multiuser communications. 

The multiuser detection system is an important aspect 
of the CDMA technology. A portion of a simplified CDMA 
detection system is indicated by the signal mixers 133 and 
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135 and the detector 137 in FIG. 1. The CDMA detection 
system plays a substantially role in terms of the system 
capacity, processing speed, interference resistance, and 
noise reduction. 

5 One significant limitation of the prior-art multiuser 

detectors in CDMA systems involves the cross-correlation 
information among users in signal detection. Many prior- 
art systems essentially neglect the cross correlation. 
This results in reduced error performance and consequently 
10 reduces the capacity of the system. This capacity loss is 
most significant when the processing gain of the system is 
high. Specifically, the cross correlation noise is not 
separated from the autocorrelated signal in many prior-art 
detectors . 

15 Another limitation of the prior-art systems is the 

slow processing speed inherent in the detectors. This is 
partially due to the complexity of the prior-art decoder 
processing techniques. For example, Verdu first disclosed 
the idea of the optimal multiuser detector in "Minimum 

20 Probability of Error for Asynchronous Gaussian Multiple 
Access Channels", IEEE Trans. Inform. Theory, Vol. IT- 
32(1), pp. 85-96, 1986. An ensemble of K received 
asynchronous signals as convolutional codes was used in 
the proposed system, where K is the number of users in the 
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system. The well-known Viterbi algorithm was used to 
decode these signals. Lupas and Verdu proved that Verdu's 
multiuser detector has a low bit error rate and is 
near-far resistant in "Near-Far Resistance of Multiuser 
Detectors in Asynchronous Channels 7 ', IEEE Trans. Comm, 
Vol.38 (4), 1990. However, Verdu's scheme exhibited two 
major problems. First, the complexity of the system grows 
exponentially with the number of transmitted signals K. 
Even for a moderate value of K, Verdu's system is already 
hard to implement. Secondly, the whole period of the PN 
sequence is assumed to be modulated by only one symbol to 
reduce the computation, which makes the trellis become 
time invariant. 

Most practical systems are usually designed with the 
whole period of the PN sequence being modulated by 
multiple data symbols. This makes the total complexity of 
Verdu's algorithm for binary phase shift keying (BPSK) 
waveforms to be 2" plus the computation of the weight of 2 
x 2 K branches in the trellis. Such complexity renders 
Verdu 1 s system even less feasible. Verdu's system also 
requires information about the received power from each 
user in order to detect the received signals. In 
addition, effects due to multipath fading were not 
considered. 
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Several suboptimal multiuser detectors have been 
introduced after Verdu's detector, aiming at reducing the 
complexity of Verdu's system. For example, Xie et al. 
described a detector using sequential decoding in 
"Multiuser Signal Detection Using Sequential Decoding", 
IEEE Trans. Commm. Vol.38 (5), pp. 578-583, 1990. This 
system has a complexity proportional to K. Lupus and 
Verdu also disclosed a linear de-correlating detector with 
complexity proportional to K. However, both of the above 
systems make the same assumption as in Verdu's system, 
i.e., the whole period of the PN sequence is modulated by 
a single symbol. 

Still another limitation of the prior-art multiuser 
detectors is the lack of capability of pipeline 
processing. The coded symbol signals are processed one at. 
a time in each channel. A new coded symbol signal is 
accepted for decoding usually after the decoding of the 
previous symbol in the same channel is completed. 

In view of the above limitations, a new type of 
multiuser detector with high capacity and fast processing 
speed is needed to improve the existing CDMA systems. In 
particular, separation of the cross correlation noise from 
the signal and reduction of the processing complexity are 
desirable . 



-7- 



WO 97/21292 



PCT/US96/I8824 



The present invention describes an optimal multiuser 
detector with a dynamical system control at a high 
processing speed. Specifically, a double-sampling front- 
end processor is implemented in the preferred embodiment 
5 in accordance with the present invention. The detector 

further includes a dynamic control module for generating a 
system control 5-parameter based on the characteristics of 
the received signals, a correlation processor for 
processing cross correlation information of the received 
) signals based on the 6-parameter, a plurality of 

approximation processors for recursively suppressing the 
cross correlation noise below a predetermined noise 
tolerance level based on the optimized 5-parameter, and a 
decision processor for generating the decoded data. The 
approximation processors form a plurality of substantially 
identical processing layers connected to each other in 
series. The correlation processor feeds information on 
both auto correlation and cross correlation of the input 
signals of multiple users to the first approximation 
processor in the series. In a preferred operation, 
received signals are processed layer by layer until the 
cross correlation noise is reduced below the predetermined 
level • 
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One aspect of the present invention is the linear 
dependence of the processing complexity with the number of 
the simultaneous users present in the system. This 
feature in combination with other unique properties of the 
5 preferred embodiment allows greater system capacity at 

faster processing speed in comparison with the prior-art 
systems . 

Another aspect is the unique capability of pipeline 
processing mechanism. Each processing channel of the 

10 preferred embodiment is pipelined to increase the 
processing efficiency and speed in addition to the 
parallel processing implementation. In a preferred 
operation mode, a layer proceeds to accept and process a 
new signal upon completion of processing the previous data 

15 while the correlation processor feeds a new control 6- 
parameter adjusted and optimized for the new signal- 

Still another aspect is the use of cross correlation 
information to substantially cancel the interference 
contributed by the cross correlation signals. This 

20 feature presents significant advantage over the prior-art 
system wherein the cross correlation signals is neglected 
as white noise. The minimized interference allows an 
increased system capacity of processing simultaneous users 
at the same and adjacent spatial locations. 
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Still another aspect is the dynamically adjusted 
system control for a high-speed convergence in the multi- 
layer processing. The detector of the preferred 
embodiment controls the decoding processing in each 
5 processing layer by feeding the layer with a dynamic 
control parameter 6. An optimized S-parameter is 
preferably adjusted and determined individually by the 
characteristics of a received signal such as the number of 
the simultaneous users, K, the processing gain, the PN 
10 sequence, and the user delay- Generation of the 5- 

parameter is a simple digital operation at a high speed in 
the present invention. 

Other aspects of the present invention include 
additional system enhancement by optimally combining the 
15 two output signals from the double-sampling front-end 

processor, digital generation of the PN sequence, and a 
mechanism for handling decoding process at very small user 
delay times relative to the signal bit duration. 

The details of the preferred embodiment of the 
20 present invention are set forth in the accompanying 

drawings and the description below. Once the details of 
the invention are known, numerous additional innovations 
and changes will become obvious to one skilled in the art. 
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BRIEF DESCR IPTION OF THE DRAWINGS 

FIG. 1 shows a block diagram of a simplified CDMA 
system. 

FIG. 2 shows the preferred embodiment of the 
5 multiuser detector of the present invention with a 
circuitry for the front-end processor. 

FIG. 3 illustrates the delay of each user signal and 
the corresponding transmitted data symbols during 
observation time from iT to (±+l)T. 
10 FIG. 4 shows the 6-adjusted Mth order detection 

module in accordance with the present invention. 
FIG. 5 shows the correlation processor. 
FIG. 6 shows the approximation processor in 
accordance with the present invention, 
is FIG. 7 shows the decision processor in accordance 

with the present invention. 

FIG. 8 illustrates the kth user's spectral-spreading 
signal, a k (t), during a observation time from iT to 
(i + l)T. 

20 FIG. 9 shows the spectral-spreading signals, a k (t}and 

a,(t) for y*" Yi during a observation time from iT to 
(i + l)T. 

FIG. 10 illustrates computation of c mn (i) for Yk>Yi- 
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FIG, 11 illustrates computation of c an (i) for Yk<Yi • 
FIG. 12 shows the process and the corresponding 
circuitry for obtaining cos(0 k -9i). 

FIG. 13 shows the magnitude of V (2k _ 1H2l _ r , , V j2k _ n (21) , 
when Yk * Yi. 

FIG. 14 shows the exact values of 6 low for p k = OdB. 
FIG. 15 shows a comparison of the exact and 
approximated values of 5 low . 

FIG. 16 shows the (3 k as a function of K/N and 5 for 
10 g k =8. 

FIG. 17 shows the P k as a function of K/N and 5 for 
g k =12. 

FIG . 18 shows the (J k as a function of K/N and 6 for 
g k =16. 

is FIG. 19 shows the (J k as a function of K/N and 6 for 

g k =20. 

FIG. 20 shows the p k as a function of 6 for K/N=0.2, 
0.4, 0,6, and 0.8 and different values of g. 

FIG. 21 shows the contour of 6 opt as a function of K/N 
20 at different g. 

FIG. 22 shows the contour of exact and approximated 
values of 6 opt as a function of K/N at different g. 
FIG. 23 shows *p opt as a function of K/N for g. 
FIG. 24 shows a comparison of 6 low and 6 opt . 
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FIG . 25 shows P .versus 6 for K/N=0.2, 0.4, 0.6, and 

0/8. 

FIG. 26 shows the results of the simulation of the 
symbol error rate probability for each user for the 
preferred embodiment of the multiuser detector. 

FIG. 27 shows the simulated average symbol error 
probability of all users versus 3K/N in units of dB. 
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DETAILED DESCRIPTION OF THE INVENT T 9H 

The preferred embodiment and examples shown 
throughout this description should be considered as 
exemplars, rather than as limitations on the present 
5 invention. In particular, the mathematical formulas, 
equations and some assumptions contained herein are 
intended to illuminate the processes involved and 
functions of various system components or are used for 
clarity of discussion and should not be construed as 
D limitations on the scope and contents of the invention. 

Many features of the present invention are described 
in "Digital Algorithms for Multiuser Detector, Decision- 
Directed Phase Tracking Loop, and A/D & D/A Converter", a 
dissertation by Lin-Lang Yang, University of Southern 
California, December, 1995, the entirety of Chapter 1 and 
corresponding appendices of which are incorporated 
herewith by reference. 

System Design of Multiuser Detector 

FIG. 2 shows the preferred embodiment 200 of the 
multiuser detector of the present invention with a 
circuitry for the front-end processor 202. The detector 
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200 is a portion of the receiving system in an 
asynchronous direct sequence spread-spectrum multiple 
access system that is also known as a CDMA system. There 
are K users operating in an additive white Gaussian noise 
(AWGN) environment. A received analog signal 201, r(t), 
at the input of the receiver is a sum of K spread-spectrum 
signals s k (t-i k ,b), wherein lsksK, and an AWGN random 
process, n(t), with two-sided power spectrum density N D /2: 



r(t) =s(t, Jb) + n(t), teR (1) 

where 

sit.b) -£ s k (t-T k ,b) = £ £b k ii)s k it-±T-x k ) (2) 

k?\ * 2-0 k*l 

io and 

s k (t) =JZvT k a k {t) cos((o c t +6^) . (3) 

In the above equations, T denotes the symbol length which 
is the same for each of the K users, s k (t) is the waveform 
transmitted by the kth user during a T-second interval, 
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a k (t) is the kth user's spectral-spreading signal which is 
normalized to have unit power, T k is the time delay of the 
kth user, 9 k is the carrier phase of the kth user, o> c is 
the common carrier frequency, w k is the power of the kth 
5 user which is assumed constant for each user during the 
observation, b k (i) is the ith transmitted data symbol of 
the kth user with b k (i)€{-l,l} for all k and i. and (L+1)T 
is the message length for each user. 

FIG. 3 shows an example of the delay of each user and 
10 the corresponding transmitted data symbols during 
observation time from iT to (i+l)T. 

A preprocessor (not shown in FIG. 2) including a 
plurality of matched filters processes the composite 
signal 201 and decomposes it into corresponding K 
is components for each user. Each component is sampled in 
the respect sampling channel in the front-end processor 
202 and converted into digital signal by analog-to-digital 
converters 210. The output digital signals from the 
processor 202 are then subsequently decoded by the 
o multiuser detection module 220 to recover the data symbols 
b k {i) from all K users embedded in the received signal 
201. 

The front-end processor 202 has K channels, one per 
user. Each channel includes a first signal mixer 204 , a 
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second signal mixer 206, and two parallel integrate-and- 
dump detectors 208 and 209. The signal mixer 204 is used 
to demodulate the received signal at a common carrier 
frequency o c generated by a local oscillator. The carrier 
phase 6 k (k=l,2,...K) for each user is determined by a 
process (which will be described herein) and used for the 
demodulation. The second mixer 206 mixes the PN sequence 
generated the a local PN sequence generator with the 
received signal in each channel for further demodulation. 
When the carrier frequency co c , the carrier phase 8 k , the 
PN sequence epoch a k (t), and the symbol time for each user 
separately are synchronized through a synchronization 
process, the received signal r(t) is despread and 
demodulated by the local PN sequence and the local 
oscillator for each user. The K demodulated outputs are 
then integrated, sampled, and converted to digital signals 
for post-end signal processing. 

According to the present invention, each channel are 
preferably sampled twice by two integration processes with 
different integration times. The first integrate-and-dump 
detector 208 for user k integrates for a duration of x k 
from time iT to time iT+T k and outputs value y 2 k-i(i)* The 
second integrate-and-dump detector 209 for the same user k 
integrates for a duration of (T-t k ) from time iT+i k to 

-17- 
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time (i+l)T and outputs value y 2k (i)- Since each channel 
generates two output signals, a total of 2K signals are 
generated for K users by the double-sampling front-end 
processor 202. One advantage of this double sampling is 
to facilitate the implementation the novel multi-order 
detection module 220 for minimizing the cross correlation 
interference and pipelined processing. 

The digital output signals y 2k (i) and y 2k -i(i) can be 
used to form a vector y(i) of 2K dimension: 



y(i) * 



y 2 (i) 
y 3 (i) 



(4) 



Accordingly, the received data symbols b k (t) embedded 
in the received signal r(t) can also be defined in form of 
a vector z (i) : 
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z{i) a 



/"A<i-1> 



<5) 



Similarly, let vector a*fi, t; represent the waveform of 
the received PN sequence corresponding to yCi>: 



a (t, i) & 



a t (t, i) 
a 2 (t, i) 



a^U, i) 
a 2jr (t, i) 



(6) 



a 2k it,i) 



a k {t) if iTs t<iT + T Jt , 
0 otherwise, 

a k (t) if ir+T^s t< (i + 1) T, 
0 otherwise, 



(7) 
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The correlation between a k (t,i) and a x (t,i) can be 



represented by c k2 (i): 



c kI (i) = a k {t,i)*a 1 (t,i) * j a k (t, i) al (t, i) dt . 



(8) 



A matrix C(i) can be formed by c^fi) to represent the 
correlation of the received PN sequence waveform: 



where a 9 (t,i) is the transpose matrix of a(t,i). 

The inventor found that ignoring the additive white 
Gaussian noise term, n(t), can simplify the description 
without affecting the f unctionability of the preferred 
embodiment. The double frequency terms can also be 
filtered by the integrate-and-dump detectors 208 and 209. 
Thus, neglect of the AWGN and elimination of the double 
frequency terms lead to the following: 



C(i) =a(t, i) •a'ft,!) 



(9) 
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Y„,<i)" J r{t)fico3{u ct + 6 k )a k it)dt 

-^^"^T (io) 

+ £ ^ J (i-l) COS ( e *- e i )C (2k- 1 )(2J-l) (i] 



2«l,2*Jk 

+ £ ^ J (i) cos < e Jt " 9 j )c (2Jt-i>m) (i) - 

2-1,2** V 



♦ £ ^ i (i-D COS ( e K- e j )C (2*M2i-l) {i) (11> 

+ £ b 1 (i)cos(e K -e i )c (2Jt)(2J) (i) . 



+ 



The above two equations can be written in vector form as 



R{i)z(i) =y(i) , < 12 > 



where matrix R(i) = [r kl (i)) is given by 
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r (2k-l) (2A-1) ' J ^ ' 



(13) 



r (2k-l) (2k) = °' 

r <2*) (2A-1) ^ 2 ' = °' 

r (2*.i)(2j-n< i ' = cos(e Jt -e i )o (2jt . 1)(2i , l) (i> , 

^-1,(2^^) = COS « 9 *- 6 i) C ( 2A.l M 2J)< i )' 

r (2A)(2i.l,« i ' = COS ( 0 ^ e 2) C (2X)(2i.l)( i )^ 



Note that r 2J (i; is independent of the data symbol index 
i, that is, r n (i) is time invariant. 

Equation (12) indicates that the received data vector 
z(i) can be decoded from the sampled digital output vector 
5 yd) by the following operation: 

z(i) =R[i)' l yli) (14) 



where R(i)' 1 is the inverse of R(i). 

Gaussian elimination is generally used to perform the 
operation of inversing a matrix. This is described, for 
example, by Strang in "Linear Algebra and Its 
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Application", 3rd edition, Harcout Brace Jovanovich, 1988. 
The matrix R(i) is a 2K x 2K matrix. The Gaussian 
elimination is known to require a number of [(2K) 3 -2K)/3 
"multiply-subtract" operations in forward elimination, and 
a number of (2K)72 "multiply-subtract" operations in 
back-substitution. Therefore, a total number of 



(2K) 3 -2K A {2K£ (15) 
3 2 



operations is required. However, this number is usually 
too large for real time processing even when 2K is 
moderate. In addition, Gaussian elimination requires at 

10 least 2K serial operations in forward elimination, and at 
least 2K serial operation in back-substitution. This 
implies that the computation time will be lower bounded 
because of the number of serial operations required. 
The inventor recognized the limitations of the 

15 conventional Gaussian elimination method and thereby 
devised a new approach to perform the matrix inverse 
operation shown in Equation (14). This new approach 
utilizes the small off-diagonal components inherent in the 
matrix R(i) to reduce the number of operations required by 

20 the Gaussian elimination. In particular, the reduced 
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number of operations in this new approach can be 
implemented using parallel processing. The new approach 
in accordance with the present invention allows real-time 
processing of z(i). 
5 The multi-order detection module 220 of FIG. 2 is 

designed to perform the operation shown in equation (14) 
to recover the data symbols represented by vector z(i) 
from the sampled digital output vector y(i). The new 
approach of obtaining R(i)' 2 is implemented. 



10 



Multi-Order Approximation Processing 

The preferred embodiment decomposes the matrix R(i) 
into two parts, an diagonal matrix D = [djj containing 
only the diagonal components of R(i) and a matrix Of 
containing only the off-diagonal components of R(i) . This 
15 can be expressed as the following: 



>* (16) 
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d <2*-l> <2*-U ' r (2*-l) <2*-l) " r 

(17) 

H - r -1-^ 

a (2*l 12*1 L <2*> i2k) T ' 



An approximation processing method is used in the 
preferred operation mode. A positive system control 
parameter 6 > 0 is introduced along with three matrices 
D 6/ 0 6 , and F 5 : 



D 6 = (1+5) D, 
0 b = 0 ( -bD, 



(18) 



5 The parameter 6 is dynamically adjusted based on the 

operation state of the CDMA system. The choice of a non- 
zero 5 is to ensure the convergence of the processing. It 
is known that matrix F 6 for 6=0 is not necessarily a 
convergent matrix. Optimization of 6 parameter will make 
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the convergent speed faster. This will become apparent 
later . 

The matrix R can be further decomposed in terms of D 6 
and F 4 : 



R= D*O f 



= (1+6) Z^I+^D^O,) j -6D 

= (l + 5)£(l-^_I + _2_ D X)) 
1*6) ^ + -±0*0,-60)) 



= ( 



= d 6 [i+d; 1 oj 

= D,(I*F 6 ) , 



(19) 



s where I is the identity matrix. The equation (12) can be 
written as 



D b (I+F 6 )z=y. (2Q) 



The preferred approximation operation can be further 
represented by a matrix H m and an estimated 
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representation of the data symbol vector z(i), Z^-M h 
both of which are defined as: 



(21) 

& H lM) y. 
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The integer, M=0, 1, 2, 3,..., is the order of 
approximation processing. 

The relation of z (i) and Z<"> is given by: 



i£o J (22) 

= (I-(-l) ml F 6 M,1 )z 

= z -(-i) H,i Fr^ 



Accordingly, the received data b^i; can also be 
represented by an estimated jb>" W • The estimated b/ w (i> 
can be determined by either z 2K m (i) or z 2K . i tw (i+D - 

A common practice in implementing CDMA systems is 
power control to ensure a substantially equal power 
distribution in the signals form K users. Implementation 
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of such power control will be assumed herein, i.e., w k 
(k=l, 2, K) for K users in Equation (3) is a constant 

denoted by w. 

The inventor recognized that the integration times 
for obtaining digital output y 2k (i) and y 2k . J (i) for the kth 
user are different and therefore different weighting 
factors should be applied in combining the output signals 
y 2k (i) and y 2k _j (i) . A combining ratio £ is therefore 
optimally obtained in implementing the approximation 
processing of Z(i). This combining ratio significantly 
affects the performance of the multiuser detector 200 in 
accordance with the present invention. The optimized 
value, also called maximal signal combining ratio, is 
achieved through a zeroth order approximation but will be 
used for all of orders in the preferred embodiment of the 
present invention. 

It will be shown that the maximal signal combining 
ratio $ is given by (T-r k )/T. Therefore, the estimated 
z 2K <M) (i) is weighted by f=(T-r k )/T, while the estimated 
Z2k-i (H) (i+l)is weighted by r k /T. Accordingly, the 

received data estimated b k (M) (i) for the kth user of the 
ith symbol can be expressed as 
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b k m U) - 



r-T, 



;(M) 

z 2k 



(23) 



where the estimation noise n M ' w is defined as function of 
FfUi-ifu™ (i)) by 



« (i) 4 -(-i^ifl^if; U)b } (i-D 

♦^fSi , , w .i J ti*i))^<i) (24) 

The estimation noise ji^'" 1 is the interference contribution 
from the cross correlation of the signals from all other 
5 (K-l) users except the kth user. Many prior-art CDMA 

systems simply neglect this estimation noise n bk m) as white 
noise in the background and consequently the resistance to 
the interference is compromised. This further adversely 
affect the maximal number of simultaneous users that the 
10 CDMA system can handle. 

The multiuser detector 200 of the preferred 
embodiment uniquely processes the cross correlation 
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information and operates to minimize the corresponding 
noise. One aspect of the present invention is to use the 
dynamically adjusted and optimized 5 parameter to ensure 
the convergence of the estimation noise n bk (H) , i.e., 
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EHn£ l) (i)) 2 ) < E{(n^(i)) 2 ) V*, i, and M, (25) 

where E{} denotes the operation of expected value. In 
addition, the 6 parameter is chosen to produce a fast 
convergence. Determination of the optimal 5 value will be 
described hereafter. The estimation noise n bk m 
contributes less variance to b k m (i) as more orders of the 
approximation processing are performed. The exact number 
of orders, M, for a particular CDMA system is determined 
by the maximal interference tolerance thereof. 

Define a new set of vectors z (M) ' wt (i)= [ z 2 w; ' wc (x) ) , x tH> , 
s w as the following: 



^ ,MMt (i) = D b (i)H< N) (i)y(i) . (26) 
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= £ (-l) m x ,m » 



(27) 



m=0 



The estimated b k m (i) can be expressed as 



b, (i) = — — z 2Jt U) + — — *2*-l U + 

(1+6) (T-Tj) ■4 2 ff (J ,,' 1 / 5 ,',l'-^^-,<i'H «•) 



(1+5)T " U*6)T ^ 

kt(1+6)V 1 



The approximation method in accordance with the 
present invention can now be revealed by the equation 
(28), which includes the following steps: 

5 1. Compute the vector z tm ' vt (i) as follows: 
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= ( £ (-D fl, (0 6 (i) D^ii) ) m j y(i) 



(29) 



m=0 



2. Determine values of b k <m (i) for the kth user as 
follows : 

if 22 r' vt (i)^ 2 , r. , r t (i + i)>o, mm b^u) -+i ; 

If z™""U)*i™i'«U+l) i0, then b^d) =-1. 



(30) 



A recursion of x ,MI and s ,M ' can be readily established 
from Equation (29) : 



(31) 



wherein 



•32- 



WO 97/21292 



rcr/usyn/iwu* 



x<°> = y, 



s< 0, = y. 



(32) 



Therefore, the estimated z <M) ' yt (i) in Mth order can be 
extracted from the (M-l)th order of s*"' 1 * according to 
Equation (29) and Equation (31). 

FIG. 4 shows a preferred implementation of the multi- 
5 order detector module 220 in the detector 200. The 

detector module 220 includes multi-order processing layers 
to implement the recursion processes in the above- 
described approximation method. The detector module 220 
includes four sections, a correlation processor 410, a 6- 

10 generator 440, an approximation processor block 420 having 
a plurality of approximation processors 422, and a 
decision processor 430. 

The 6-generator 440 monitors the operation status of 
the CDMA system and produces an optimized 6 for each pair 

15 of the input signals 412 and 414 to the zeroth order 

approximation processor (in fact, they are identical to 
each other) . The optimized 6 to the correlation processor 
410 automatically changes with operation status of the 
CDMA system. This dynamic mechanism is to keep the 

20 processing efficiency and speed at optimal. 
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The correlation processor 410 processes cross 
correlation information based on the PN sequence a(t,i) 
402, relative phase difference between users 404, i.e., 
cos (6 k -&j) , and the dynamically adjusted 6 value 442. 
5 FIG. 5 further illustrates the main tasks performed by the 
correlation processor 410 including computing the matrix 
R(i), and computing the matrix D 6 . The information of the 
autocorrelation is included in the output matrix D 6 and 
the information of the cross correlation of different 
10 users is included in the output matrix 0 6 . Both output 
signals D 3 and 0 6 are 6-dependent and thereby are 
dynamically adjusted to the status change of the CDMA 
system. 

The output signals from the correlation processor 410 
15 are fed to the approximation processor block 420 to 

generate an output S tH} 424 with a minimized interference 
due to cross correlation of different users. A plurality 
of substantially identical approximation processors 422 
are connected in series relative to each other to form 
20 multiple processing stages. As described previously, the 
5 is so chosen to ensure fast convergence of the multi- 
order processing so that the estimation noise n tK m> 
decreases as M increases. FIG. 6 further illustrates the 
operation of the approximation processor 422. 
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One unique feature of the approximation processor 
block 420 is the capability of continuous processing in a 
pipeline fashion. Each processor 422 at jth stage accepts 
and starts processing a new set of data from (j-l)th stage 
5 upon completion of feed previous set of data from (j-l)th 
stage to the next stage (j+1) . Few prior-art systems 
implement such pipeline operation. The pipeline operation 
allows significant enhancement in processing efficiency 
and speed, thereby making the real-time implementation of 
10 multiuser detector feasible. 

The approximation processor block 420 feeds the 
resulting signal s (M) to the decision processor 430 if the 
noise level from the cross correlation is reduced below a 
pre-set tolerance level. The operation of the decision 
15 processor 430 is illustrated in FIG. 7. 

The above-described processors and the operations 
thereof can be implemented with the existing electronic 
circuitry and VLSI technology. 



Correlation Processing and Correlation Matrix R 



20 Signal decoding and detection by correlation 

operation is a known and widely used processing technique. 
The complexity of the correlation process usually 

-35- 



WO 97/21292 PCT/US96/18824 



increases significantly with the number of simultaneous 
users present in a CDMA system. The necessary processing 
or computation often makes real-time processing 
impractical with the current digital technology. The 
preferred embodiment uses a new and unobvious approach to 
reduce the complexity of the such correlation processing 
by processing the auto correlation and cross correlation 
in a unique way. This is partially shown in the 
correlation matrix R(i) that organizes the correlation 
information for further processing. 

The correlation matrix R(i) is defined by Equations 
(12) and (13). Information about the matrix C(i) and 
cosfe^dj) are needed in order to obtain the matrix R(i) 
at the receiver. 

In CDMA technology, a symbol bit is spectrally spread 
by the spectral-spreading signal waveform of a square wave 
Pit). A single pulse of the PN waveform p(t) is called a 
chip. The chip duration T c determines the symbol 
processing gain W: 



Processing gain N = — 

C 
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The spectral spreading signal for the kth user (normalized 
to have unit power), a k (t), over a T-second interval can 
be expressed as the following: 

a k (t) =£±u k {n)p{t-nT c -Y k ) (34) 

n*0 Jf 



where Mn)is the nth transmitted chip of the kth user 
5 with u k (n)€{-l,l) and y k is the chip delay of the kth user 
where \>,=x k modT c . The waveform a* ft; is illustrated in 
FIG. 8. 

The computation of the correlation matrix C(i) is 
preferably obtained digitally in accordance with the 

io present invention. This is due to the special structure 
of the above spectral-spreading signal a k (t). A digital 
circuitry suitable for generating such function is well- 
known in the art. Digital generation provides higher 
stability than the analog counterparts at lower cost. 

15 Let the chip sequence, sh k (i), of length N+l for the 

kth user between time iT and (i+l)T be defined by 

sh x (i) = {u t (iN-l), u k iiN), (iN+1) , u, ( (i+1) N-l ) ) . (35) 
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In addition, let sh k * T (i) denote the right-shif ted version 
of sh k (i) , that is, 

sh k (i) -*1 

(36) 

iu k (iN-2) f u k (iN-l), u k (iN) , u k (ii + )N-2)\ 

where ->j denotes right-shift by j chips. If the chip 
delays Yk an d Yi are both the same as shown in FIG. 9, 
5 then it is straightforward that c (2k , }> (21 - 1} (i) t c (2k , v (21) (i) , 
c (2k) (2i-i) d) r and c (2k)(2i)d) are functions of sh k (i) esh x (i) , 
T k and T lf where e> denotes the exclusive or operator. The 
calculation of sh k (i) &sh } (i) can be implemented in VLSI 
technology. This is well known to the art. 

10 Furthermore, most of the computation involved in the 

above derivation is merely the sum of random +1 and -1 
plus the fractional chip correlation at the edge. The sum 
of random +1 and -1 only takes a very small amount of 
computation time in digital signal processing. This is 

15 one reason the preferred embodiment uses digital 
implementation for generating R. 

In the case that Yk >Yh the preferred embodiment 
first operates to move a k (t) leftwards in time domain by 
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an amount of (Yk ~Yi) to achieve chip alignment. This is 
to calculate the correlation values digitally. In this 
case, the correlation thus obtained is represented by 



with m = 2k-l, 2k and n - 21-1, 21. The variable c mn LT (i) 
is a function of sh>, (i) ©sh A (i) . 

Secondly, the preferred embodiment first operates to 
shift a k (t) rightwards in time by an amount of T c -(y* ~Yi) 
and the corresponding correlation thus obtained is 
represented by 



<CU), (38) 

io in this case, c nn RT (i) is a function of sh/ r (i) Qshj (i) . It 
is observed that the exact correlation value is given by 



c„li) -a w ei T (i) + (l-a w )ci r (i) (39) 



where a kl is given by 
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(40) 



FIG. 10 illustrates the derivation of c an (i) for y k >Yi . 

Similarly, when Yk < Yi/ a 2 (t) is shifted first 
leftwards and then rightwards instead of moving a k (t). 
The corresponding correlation values thus obtained are 
denoted by c mn LT (i) which is a function of sh, (i) esh k (i) and 
Cm* T (±) which is a function of sh/ 7 (i) &sh k (i) , 
respectively. Accordingly, the correlation value is given 
by 

^n c ^d) + (l-a w )c£f(i) (41) 




(42) 



This process is illustrated in FIG . 11. 

It is concluded from Equations (40) -(42) that the 
correlation value can be generalized as 
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Cli) = a kJ c£(i) + (l-a kl ) c£{i) (43) 



with m = 2k-l, 2k and n = 21-1, 21, where a kl has value 
given by 



T T 

c c 



(44) 



In addition, the ones with larger chip delay shall be 
moved rightwards and leftwards. The cornerstone of the 
5 proposed multiuser detector is that the derivation of 

c mn LT (i) and c mn RT (i) mostly are the sum of random +1 and -1. 
This particular feature makes the characteristics of the 
matrix R(i) readily obtainable digitally. Furthermore, a 
method of obtaining an approximation to the inverse of 
10 R(i) can be determined based on the characteristics of 
R(i) . 

FIG . 12 illustrates the circuitry and method for 
generation of cos (0 k -6 x ) that are implemented in the 
preferred embodiment. A reference signal at the local 
15 oscillator frequency co c with a common reference phase 9 0 

(an arbitrary real number) is generated for all K users in 
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the CDMA system. It is again assumed that the integration 
process will filter out all the terms with double 
frequency. 

For the purpose of more concise notation, the system 
5 time index will be omitted mostly thereafter. Moreover, 
only time invariant symbol delays are considered. 

Define Vmn as the time overlap between a n (t,i) and 
a m (t,i) in unit of chips. FIG. 13 illustrates examples of 

V (2k-1) (21-1) t V (2k-l)(21)' V (2k) (21-1) anC * V (2k) (21) • 

1° From Equation (13), it is observed that with time 

invariant symbol delay, the components r (2k . u (2k . 1} , r (2ki(2k) , 
r (2k-D(2k)t and ?(2k)(2k-i) f°r k=l to K, are all constants with 
values of 



<2A-1) (2*-l) 
r (2*> (2k) 

T (2k) (2*-l> 



r 

T V 

2_ * _ "(2*1(2*) 



{2k-l) [2k-l) 
N 



N 



(45) 



0, 
0 . 



Whereas, Equation (13) and Equation (43) indicate 

that r ( 2k-i) (21-1) is a random variable with a value given by 



.42- 



WO 97/21292 PCT/US96/1 8824 



r <2*-lH21-l) = COS (6 A -6) C (2/c . 1)(2i . X) 

I LT RT \ <46) 

= COS (9^-6) \<* kI C {2 k-l\ (2J-1) + < 1 ~ 0C A J E ) C 2*-U <2J-l)/« 

In general, the probability distribution of 9 k and a kl , and 
the relationship among random variables 6 k , 0), c m -iu2i-i)* 
<* kl , c LT (2k . 1)(2J ^, c* T (2k . 1)m . 1} can be summarized as follows: 

1. The random variable 6 k and B x are independent and 
5 identically distributed with uniform distribution between 

-n and n, and they are independent of c (2k . u (2 i-i)» 

2. The random variable a kl is uniformly distributed 
between 0 and 1, and is independent of c LT (2k . }) (2i . u and 

RT 

c <2k-U (21-1) • 

10 The moments of cos (e k -e t ) and ot w can easily be 

calculated with the above relation: 



EtcosO^-ej)} = 0, 

E(cos 2 (e t -e.)} = , 

2 

ElcosMe.-ep) = 1, 

I, 



(47) 
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Assuming the chips in a PN sequence is independent 
identical distribution ("i.i.d") with equal probability to 
be +1 or -I, then the moments of c (2k _ n (21 _ J} LT and c HT (2k . ln2I ^ 
can be derived as shown in the following: 



E( c i2k-l) (21-1)} ~ E ^ C 2k-l) (21-1)^ = 0/ 
£ ' ( C ( L 2*-1> {22-D ) ^ s E ^ C 2*-l) (21-1) ) ^ ~~2 V <2k-l) {22-D ' 



£ '' C (2ic-l)(2J-l)) 4 ' = t C 2Jc-l) <2 J-l ) ) ^ 



N 



1 2 

7 < 3 V 2k-l) i21-l)~ 2V {2 k-\) (2J-1) > • 



and 



^ lc (2lr-l) (21-1) c (2*-l) (21-1)' " u / 



(48) 



1 I C (2*-l) (21-1) * i C (2Mj(2J-l)l ' - ' 

' ( C (2*-l) (21-1) ) 1 ( C (*2*-l) (21-1) ) ^ = c (2k-l) (2J-1) ) 3 ( C (2A-1) (21-1) ) 1 



(49) 



= 0, 

Therefore, the moments of r m -m2i-n ar *d r lk are given by 
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and 



2 1 
£ ' r (2*-l)(2i-l)' = ~^~2 V l2k-\) i21-l) ' (50) 



3AT 



I* V 1 [ 54 t/ 2 3 V I- 

£^(2*-lM2i-l)' s ^7 \-Jq V (2*-l)|2i-l) Yq v <2k-l»(21-l)/' 



Matrix F 6 represents information of the off-diagonal 
elements of the correlation matrix R(i) . F 6 is 6- 
dependent and allows optimal control of the correlation 
5 processing in minimizing the interference- Each 
components in F 6 and F 6 2 can be expressed as 
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- (1) 6 

'A" ■-— ^f* 1 * £ -^-). (52) 
(1+6) 2 I "lr w rjl / 

The first and second moment of f kk m , f k i <l: , f^' 2 ', and f kl t21 
can be computed as follows: 
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tit") - - -r—r 
1+6 

Elftf 1 ) = 0, 



kk (1 ♦ 6) 2 V 3 v u v j/ 

E{(f = $!_, 

kk (1 ♦ 5) 2 

. _^ ((*» . f £ i £)' 

kk (1 +6) 4 H V kk j.i.j.k 3 V w / 

_L ? filial) 



The following relation is observed in obtaining the above 
results : 

Zirfri) - E(r£te(r£j. (54) 
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Determination of Maximal Signal-Combining Ratio 

The maximal signal combining algorithm dictates the 
performance of the proposed system. Due to the complexity 
that is involved in deriving the combining ratio, the 



zeroth order multiuser detector. The same combining ratio 
will be used for other orders of multiuser detector. 
Although such weighting ratio does not guarantee the 
maximal estimation signal-to-noise ratio for M *1, it 
avoids the totally intractable computation that can be 
involved. Furthermore, as M increases, the estimation 
noise will be dominant by the adjustment 6 instead of the 
combining ratio therefore, it make sense to use the 
zeroth order combining ratio for all of the orders. 
Suppose the estimation b k iH) (i) is generated by 



5 



maximal signal combining ratio is calculated only for the 




(55) 



where the estimation noise n t>k ' m (i) is given by 
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♦ £(tfffi , {i 1 ,u)*(i-t)f 1 si , 1 «..,(i*i))«' J (ii < 5S > 

3 = 1 X 

* H-^«, 2j , (i + l)b j( i + l)). 



Then the zeroth order multiuser detector will have 
estimation signal-to-noise ratio, SNR<°>- " (i) , given by 

E(jb5(i)> l 
SNK«»'*U> • — -"J -137' < 57) 



(2k) (2j) 



(i) ) 2 ) 



♦ 2U1 " ?)Elf£!M2j) (i))E(f { , 2 1 k , .xH2j-ii < i + 1 >> 
+ E{((l-Of/2k-U( 2j -l)(i + 1 )) 2} 

♦ E(((l-l;)«,2k.in2j) (i + D) 2 >) 

_L_( 6 ' + I*±( < 2 ^ ^ 



(1 



(2k) (2*) 

N 



{2k) (2k) 
N 



(58) 
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Maximizing SNR (0) ' k (i) is equivalent to minimizing p Ui) 
which in turn is equivalent to minimizing the following 
term: 

V i2k) (2k) ^ _ V (Zk) Uk) (59) 
N N 

which reaches its minimum value at £ max : 

This is the same signal combining ratio that is used in 
Equations (23) and (24), i.e., (T-r k )/T. 

Although a similar procedure can be done for the 
second order multiuser detector to find the maximal signal 
combining ratio, it turns out that the combining ratio is 
a very complicated function of V kk for i = 1, 2, 3,..., 
and 2K. Likewise, it is infeasible to calculate the 
precise v (M) for M *2 which implies that the exact maximal 
combining ratio is very difficult to obtain for M *2. 
Therefore, the maximal combining ratio of the zeroth order 
detector will be used for other orders as well. 
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Moreover, as M increases, the multiuser detector 200 
relies on the adjustment 6 to reduce the impact from 
n bk (H) (i) instead of the combining ratio 

Estimation Noise Variance o tM,,k and Optimization of 6 

As discussed previously, the estimated b k m (i) and 
the estimation noise n bk ' M > (i) are related by b„ ,M > (i)=b k (i) + 
n bk <H) (i) given by 



J -1 -t 



a <M),k is defined as the estimation noise variance for 
user k that is generated by the Mth order multiuser 
10 detector: 
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,M >'* £E{(n b ( t M, (i)) 2 } 

=e{(-(-1)- 1 '£ (J^ (D^d-i, 

♦(-^ f ffw/u, (i) ^fK„2 H , (i) 

4 f K>^ li+llij ; (i+1) ))} 

=£ ( f .sr,v„.„ (i) ) 2 } * E Jill* f «£» „ {i) ) 2 } 



where cr' w -* denotes the estimation noise variance 
contributed by user j to user k in the Mth order 
multiuser detector and is given by 



. e{( I^f <i.) 2 } . e {( l^ifffi",, u,) 2 } 

* 2 T "^ Elf '"'«"' U1 ,E,f «?"> <»-., (i * i) ) 

♦«J^'.S?i«..,ii*i))'} 

^{(^f,K„ 2j> (i*l)) 2 }. 



The order of the expectation and the summation in Equation 
(62) is exchangeable, because the received data of 
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different users are independent relative to each other. 
Also, it is assumed that the received data are equally 
likely to be +1 or -1, that the random variable f kl m (i) is 
independent to the received data bjU), and that the 
random matrix F 6 ,M1 (i) is independent of the random matrix 

F 6 ,MI (1) when i * 1 . 

Since the first and second moments of F s and F 4 2 are 
available in Equation (53), a^ 0 ''" and o 5 (1 '' k can be derived 
as the following: 



(1),* 



(1 


+ 5) 2 


1 


1 


3N 


(1 ♦ 




1 


(1 






1 


(1 


♦ 5) 4 


1 

+ — 


U 


N 





(64) 



io wherein 
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Special attention should be given to the 
user-independent second moment of the estimation noise, 
a/ 0 ''" and a/ J ''*, when 6=0. This particular structure 
happens when the proposed maximal signal combining 
5 algorithm is used. 

The optimal value of 6 for order M is the one that 
achieves the minimum estimation noise variance from all 
users, i.e., 
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6 opt 



= U :min£ o< w >'4- < 66 > 

I 6 *»1 J 



Equation (64) indicates that the complexity of the 
estimation noise variance as a function of 6 increases as 
the order increases. Therefore, it is very difficult to 
find an optimal value of 5 for all orders. 
5 The multiuser detector 200 of preferred embodiment 

instead searches for a 6 which assures the fast 
convergence of the estimation noise. Moreover, the 
optimal value of 5 should not be order-dependent or 
symbol-delay-dependent in order to reduce the 
10 implementation complexity. 

Using Equation (62), it is found that the upper bound 
to the estimation noise variance is given by 
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^{(-^£S--V, l2J -„(in)) 2 } 
^(^fS.V.^.iin)) 2 } 



wherein E(f kl iH) )= 0 for k * 1. 

When M is a moderately large number, the first term 
in Equation (67) will be much smaller than ~he second term 
when the estimation noise variance is a decreasing 
5 function of 6, i.e., 



r- 1 



« O 



j*l i-1 J 



(68) 
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Therefore, the upper bound of the estimation noise 
variance o (M, * k can be approximated by 



a m.k . o (0) # *g gdf™}. (69) 
u upper £J Hi iJ 



Similarly, we can also derive the upper bound to the 
estimation noise variance a 4 "* 1 ' 0 *: 



a W+l),k m 0 {l),*g £ (70) 
upper JJ 



It is desirable to find a 6 value so that 



= 16 : min ), 



6 = 15 : min ), (71) 



° pt * E^ 0 o ,M, ' A 



which will allow the optimal decrease of estimation noise 
between two consecutive processing stages in the 
approximation processor block 420 shown in FIG. 4. 



•57- 



« 



WO 97/21292 



PCIYUS96/18824 



Since the close form of o (H> ' k is very difficult to 
find, the 6-generator 440 shown in FIG. 4 is thus 
configured to obtain the optimal 6 opt by the following 
criterion instead, 



5 t = is : min Z *«°<«>" ) 

^k^upper /p ^ % 



6 : min- 



6 r* n <0),* 



5 The apability of serving a large number of 

simultaneous users is one of the advantages of a CDMA 
system. Therefore, it is reasonable to assume that the 
number of users, K, is much larger than 1, i.e., K » 1. 
This assumption is used in the description herebelow. 
1 <> The ratio of the estimation noise variance of the 

first order to that of the zeroth order for the kth user 
can be found as 



3* * 



r(0) t k 



(73) 



WO 97/21292 



PCT/US96/18824 



} <D.k 



K lll.k 

X* o i0) ' k 



(6 2 + — — ) 2 + — ( — ) + — - 5 2 — 
3 AT 9 N 3 N 

(6 2 + 1JZ) (1 + 6) 2 
J N 



(74) 



wherein 



J s i f g[ V M-»l2 k -l) V (2j-l)(2;j-l)) (?5) 

It is desirable that (J k is less than 1 for all k in 
order to assure convergence of the preferred multiuser 
detector 200 in FIG. 2. 
5 For the general case, 6 shall satisfy: 
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5 3 +( i- JL-1^,6 2 + 1 + l * ( *>' > 0. 
2 2N 6 AT 3N6W3W 



The lower limit for 6 is approximately given by 

K 1 if — ,1 . _ _ k 

6, - — + gJc - . (77) 

IOW 2N 6 A7 2 

FIG. 14 shows the exact value of 5 low for g k =8, 12, 
16, and 20. The exact value and the approximated value of 
6 low are compared in FIG. 15 for g>12 and 20. 
5 FIGs 16-19 shows 3 k (in dB) against K/N and 6 for 

g k =8, 12, 16, and 20, respectively. It can be seen that 
the area above 0 dB increases with g" k , especially for 
small 6 and large K/N. 

FIG. 20 shows the values of (J k (in dB) against 6 for 
10 K/N=0.2, 0.4, 0.6 f and 0.8. 

It should be noted that g* k be bounded. Otherwise, 3k 
will easily exceed OdB for large value of "g k . This can be 
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understood by the result that the expectation value 
E{g(x,y)} is large when x is close to zero and g k will be 
close to E{g(x,y)} for x= V uk ., (2k . 1) /N when K is large. 
Therefore, one of the alternatives for avoiding a large g„ 
is to discard the information carried by very small V n . 

For example, the information carried by V u < eN 
where e«l . In this case, the 1th row and the 1th column 
of matrices O s and Dj' 1 are set to zeros, and the 1th 
component of vector y(i) is set to zero. Moreover, the 
combining ratio, i, in Equation (55) will be zero if 1 is 
an even number and 1 if 1 is an odd number. 

Therefore, in a preferred operation of the multiuser 
detector 200 in FIG. 2, the output signals, y 2k .,(i)/ from 
the integrate-and-dump detector 209 will be set to zero 
15 when the delay x k of the kth user is smaller than a 

threshold value, T min , since the filtering effect of the 
integration sampling process is substantially eliminated 
when i k <T mln . One example of setting the threshold value 



10 



*«i» x -' (78) 

nun £j r 



20 where T is the symbol bit duration and N is the processing 
gain . 
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The optimal 6 value produces the minimum average 
value of 3k tor all K users: 



3 opt i min3 



= min A £ P A 

6 A **1 



3 W 9 W 3 AT 
mm . 

• (6 2 ♦ 1£) (1 ♦ 5) 2 

3 N 



9*±ig k . (80) 



which results in the following equation for 5 opt : 



J opt 



- (1 + _ g ,6 opt + -(-)5 0 . 



opt 



K 2 



2 K 3 
9 W 



(81) 



FIG. 21 shows 6 opc vs. K/N for different values of g 
It can be shown that 6 opt can be approximated by the 
following formula : 
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5 opt - | (1 * jg) (82) 

The exact value and the approximated value of 6 opc are 
compared in FIG . 22, which shows a small deviation. 

The preferred embodiment of the present invention 
uses the approximated value of 6 opt of Equation (82) as a 

5 criterion for choosing an optimal the approximated value 
of *6 opt . Parameter g in Equation (82) usually does not 
change significantly. Hence, the optimal 6 opC is 
essentially a function of number of users, K, and the 
processing gain, N. 

10 FIG. 23 shows the average ]J opt in dB as a function of 

K/N for different values of g. It is shown that (i opc can 
be as small as -5.5 dB for K/N=0.15. When K/N=l.l, *3 opt 
still has a value of -ldB. 

FIG. 2 4 compares 6 opt and 6 low in terms of K/N for 

15 different values of g. It is shown that 5 opt satisfies the 
convergence requirement for g k *16. 

FIG. 25 further shows that 6 opt indeed achieves the 
optimized value 3 opt f° r K/N=0.2, 0.4, 0.6, and 0.8. 
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Simulation Results for Preferred Multiuser Detector 



Several operation scenarios of the preferred 
multiuser detector 200 were stimulated to demonstrate the 
some aspects of the present invention. 

It is assumed that each users 1 PN sequence is 
generated by using the following generator function. 



g(x) =x 42 +x 35 +x 33 +x 31 +x 27 +x 26 

+ x 25 +x 22 +x 21 +x 19 +x 18 +x 17 +x 16 (83) 

+x 10 +x 7 +x 6 +x 5 +x 4 +x 3 +x 2 +x + l 



The symbol delay, T k , the received carrier phase, 6 k/ and 
the transmitted data, b k (i) f for all users are randomly 
generated. The input signal y(i) to the multi-order 
detector module 220 of FIG. 2 is obtained by using 
R(i)z(i)=y(i) . 

A symbol data rate is subsequently generated by 
comparing the estimated data, b k {M) (i) with the transmitted 
data, b k (i) . FIG. 26 shows the symbol error rate of each 
user for the 6-adjusted Mth-order multiuser detector 200. 
Comparison are made for M^O, 2, and 4. Each data point is 
generated with 10* simulated symbols. In this simulation, 
the number of users in the system, K, is set to 40, and 
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10 



the symbol processing gain, N, is 42. Moreover, the e is 
set to 1/42=0.0238, i.e., any information carried in any , 
interval shorter than 1 chip will be discarded. It can be 
shown that the variable g has a value of 6.9 for 
6=1/42=0.0238. The optimal value of 5 is approximately 
given by 6 opt = (K/N) ( l+g/3) = 3.143 for g"=6.9. 

FIG. 27 illustrates the average symbol error rate of 
all users against 3N/K in dB for M = 0, 2, 4, and 6. The 
parameters, except K and 6 opt , are set to the same values 
as in the first simulation shown in FIG. 26. Again, each 
point is generated with 10* simulated symbols. 

FIG. 27 indicates that the system capacity increases 
as M increases. The system capacity is usually defined as 
the number of users that can be accommodated in a system 
15 with the average symbol error rate below some threshold 
value. For example, if the performance target is set to 
keep the average symbol error rate below 10" 2 , then FIG. 
27 shows that there is a capacity gain of about 2.5 dB for 
M = 2 as compared to M = 0. 

The above described multiuser detector can also work 
for a system with inclusion of the background noise such 
as additive white Gaussian noise. 

In addition, the present invention can also used to 
improve incoherent system. Many terrestrial wireless 



20 
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applications use non-coherent detection because each user 
has independent random phase due to Rayleigh fading 
channel. The coherent detection can easily be used in 
forward link (base station to mobile station) detection in 
an IS-95 based CDMA cellular system because the existence 
of the pilot signal. The multiuser detector of the 
present invention can be used to cancel interference 
caused by signal from adjacent cells, especially when the 
mobile station is in a hand-off region. 

A number of embodiments of the present invention have 
been described. Nevertheless, it will be understood that 
various modifications may be made without departing from 
the spirit and scope of the invention. Accordingly, it is 
to be understood that the invention is not to be limited 
by the specific illustrated embodiment, but only by the 
scope of the appended claims. 
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CLAIMS 



What is claimed is: 

1. An information processing system for detecting 

encoded signals of a plurality of users in a CDMA 
system, comprising : 

a sampling processor, having a demodulating 
element and an integrating element for sampling 
said encoded signals, producing a plurality of 
digital signals indicative of said encoded 
signals/ and 

a detector, receiving signals from said sampling 
processor, having a correlator for performing a 
correlation of said encoded signals to obtain 
indicia of interference of said signals from 
different users and a processor for processing 
said indicia of interference in a predetermined 
way to substantially remove said interference 
from said detector. 
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A system as in claim 1, wherein said integrating 
element has a circuitry for sampling each of said 
encoded signals with at least two integration times. 



A system as in claim 1, wherein said processor of 
said detector includes a plurality of processing 
modules connected in series, operating to 
sequentially process said encoded signals based on an 
approximation process, each of said processing 
modules responding to a system status parameter in a 
dynamic manner. 

A method for minimizing interference of cross 
correlation of received signals from different users 
in a CDMA system, comprising: 

sampling said received signals to obtain 
multiple representations of said signals based 
on coding information thereof; 

processing said signals in said representations 
with a first correlation operation; 
processing said signals with a second 
correlation operation subsequent to said first 
correlation operation, said second correlation 
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operation operating according to a operating 
condition of said CDMA system; and 
performing a plurality of recursive operations 
to minimize said interference. 

5. A signal improving device for a CDMA system, 
comprising : 

a signal processing element, obtaining a CDMA signal, 
and processing said signal to determine a portion of said 
signal which represents interference; 

a filter, removing at least a part of said portion of 
said signal that represents interference, to produce an 
interference-filtered signal; and 

a recursion device, receiving said interference 
filtered signal, determining an amount of interference in 
said interference filtered signal, and sending said 
interference filtered signal through said signal . 
processing element and said filter if an amount of 
interference in said interference filtered signal meets a 
specified criterion. 
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The J-adjusted Mth order multiuser detector. 



WO 97/21292 ^ PCT/US96/18824 



a(U) 



-\ 
-/ 



co$(9 k -8 1 ) 







D 5 N 


(i) compute matrix R, R = D+O f 






(ii)D s = (l+5)D, 


O s = O r 6D. 












JL MM2_ 





CCB(Correlation Computation Block) 

The correlation computation block. 



Mio 



Figure 5 



WO 97/21292 PCT/US96/I8824 

6/27 



6.02.- 



-\ 
-/ 



s <m> (0 



-\ 
-/ 



x (M + v (i) = o s p-ix lM >(i) 



s (M+l) (i) = s <M) (i) + M yu+l ^W f() 




6*12- 




■072. 



/I 



APB(approximation processor block) 



Tie Approximation processor block. 



Figure 6 



7/27 



PCTYUS96/18824 



if(iSi)+s ( ^. 1 )(i+D>0)( 

decide b$$ = +l; 
} 

else { 

A(M) 

decide b k (i) = -1; 
} 



£(M) 



DB(Itecision B/ocfy 

The decision block. ^\30 



Figure 7 



WO 97/21292 PCT/US96/18824 

8/27 



IT 



1 

Tt 



a 2i _,('.*) 



-1) 



u A (fAT+l) 



a 2Jk (f,0 



MO 



1) 



The ah user's spectral-spreading signal, «»(*), during the observation 
time »T to (t + 1)T. 



Figure 8 



WO 97/21292 PCI7US96/18824 

9/27 



iT 



«,»(«. A 

>• 



(•+1)7- 



The spectral-spreading signals a k {t) and aj(0 when 7k = H during 
the observation time iT to (i + l)T. 



Figure 9 



WO 97/21292 PCT/US96/18824 

10/27 



IT 



, («♦ 1)1- 



move 
rightwards 
by 



ST 



<•) 



* « 21-1 « »' 




Computation of (^(0 when 1 k > 



Figure 10 



WO 97/21292 PCT/US96/18824 

11/27 



Lii leftwards 




„ . or it 

a m II- I " 2 ' 



Computation of Cm„(t) when 7k < T»- 



Figure 11 



WO 97/21292 



12/27 



PCT/US96/18824 



rr 



4cos(a» c r*0,) 



0— 



integrated dump 




yCOS (CO c /+e 0 ) 



co$(0 o -ej) 



integrated dump 



^jin(o) c r + e 0 ) 



$in(8 0 -e,) 



0— 



integrate k dump 



"7V 



^co$(u) c f+e 0 ) 



cos(e 0 -e 2 ) 



0— 



integrated dump 



^$in(a) c f + e 0 ) 



$in(e 0 -e 2 ) 



0— 



Integrate d dump 



:£cos(ay + 8 0 ) 



co$(e 0 -ep 



iT 



0— 



integrated dump 




sin (o 0 -ep 



I 

o 

o 



© 
i 

o 



s 



cosce^-e,) 



The derivation of cos(0* - 



Figure 12 



_ mmmM PCT/US96/18824 
WO 97/21292 

13/27 



IT , 



V (tt-l)<2/-l) 



IT , 



fl 21 (U) 



V (2i) (20 



f W (2J-1) 



The magnitude of V ( 2*-i)(2/-i)> VJafc-i)P0» V)(^-D and V WW when 



Figure 13 



WO 97/21292 



14/27 



PCT/US96/18824 




Figure 14 




Figure 15 



PCT/US96/18824 

WO 97/21292 

16/27 




Figure 16 



r 



PCT/US96/18824 

WO 97/21292 

17/27 




Figure 17 



WO 97/21292 



18/27 



PCT/US96/18824 




Figure 18 



WO 97/21292 



19/27 



PCT/US96/18824 




Figure 19 



WO 97/21292 



20/27 



PCT/US96/18824 



K/N=0.2 



K/N=0.4 



fik 
(dB) 



ft 
(dB) 



-1 




-3 




The A in dB against * for K//V = 0.2,0.4,0.6 and 0.8. The four 
L „ are for 9l = 4 (solid line), & = 8 (- - -), * = 12 (-), 



curves from the bottom up are for g k = 4 (solid line), 
and gu = 16 (--), respectively. 



Figure 20 



WO 97/21292 



21/27 



PCT/US96/18824 




Figure 21 



PCT/US96/I8824 



WO 97/21292 

22/27 



5 op t 2 




1 1.1 



The contour for the exact value and the approximated value of 
The four curves from the bottom up are the exact contour with g = 6.9 (soud 
line), the approximated contour with § = 6.9 (...), the exact contour with g = 7.9 
(. - -) and the approximated contour with g = 73 (-.-), respectively. 



Figure 22 



PCT/US96/18824 

WO 97/21292 

23/27 




Figure 23 



PCT/US96/18824 

WO 97/21292 

24/27 




Figure 24 



WO 97/21292 



25/27 



PCT/US96/18824 



K/N=0.2 



K/N=0.4 






1i 








(dB) 


0 


VN : 
\' : .\ : 


,...,«...*•• 


». . • ' 




-1 


....\-A.i 

%. 


■>»••••••••• 













2 




(dB) 







K/N=0.8 


3 






2 






1 












0 




" • • • 








-1 




Ss^/^. . > 



The j5 in dB against 5 for /T/N = 0.2,0.4,0.6, and 0.8. Thefeur 
curves from the bottom up are for g = 6.4 (solid line), g = 6.9 (- - -). 9 - 7-4 
(...), and g = 7.9 (-.-), respectively. 



Figure 25 



* 

PCT/US96/18824 

WO 97/21292 

26/27 



10 

...::::::L::: ! :.:::::i::H:::MM=S^H>::n/!||n!l!!!l 

:.-v**::!:::::::::i:;:::::sa..;.i....."- 
% *Vc *** V"" 



o 



10 



10' 




B ...»!HjiH!!i!i!ii!iHi!!!iiaai*i!!Hi!ifaHi!!!!iiij!!!H!ij!!!!!!H 

; ; f 



"" :::: v";!;;i;!i!!i:i:::::::::-::::::::::-::::::::::- 



J W 15 20 25 30 35 40 

which user 

Symbol error probability for each user. 



Figure 26 



PCT/US96/I8824 

WO 97/21292 

27/27 



10 



10" 














3 

I 


% ! 



Eiii-:::::.^ 

••••*a : 

\ . 


2 i 


K'\ 


»— 
O 




o 






J 



TTTTTT 



r- 

0»'**y ** 



•7*0;* 



- • ™.r. :::::: - 

:;:::::::r^mM 17.1v, '.v.v.'.iv.'.'.v.'. ........... 

* r * * * ' ! \ \ * I ! \ I * ft vl! ! ; ! ! ! \ 1 1 ! I \ ! ! ! " 

... ~q ..•...»•••••.•••"• •..».... * ^ 

Ni •> 5 ::1:::::::.. 

\ iL ::::: L ;:::::: J 

'"t^! I"""*'* ••••••••• > 

.V... s« r ; r • 

, : 

....\.; ^*-v^ ; 

....A < V'-:- ••• w * 

V : i 
■: : 

i I i 



5 



1 — 

fin dB 



10 



Average symbol error probability of all users versus ^ in dB. 
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